Quantum Smoluchowski equation: A systematic study 
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The strong friction regime at low temperatures is analyzed systematically starting from the for- 
mally exact path integral expression for the reduced dynamics. This quantum Smoluchowski regime 
allows for a type of semiclassical treatment in the inverse friction strength so that higher order 
quantum corrections to the original quantum Smoluchowski equation [PRL 87, 086802 (2001), PRL 
101, 11903 (2008)] can be derived. Drift and diffusion coefficients are determined by the equilibrium 
distribution in position and are directly related to the corresponding action of extremal paths and 
fluctuations around them. It is shown that the inclusion of higher order corrections reproduces the 
quantum enhancement above crossover for the decay rate out of a metastable well exactly. 
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I. INTRODUCTION 

Real systems interact with surrounding reservoirs 
which typically contain a macroscopic number of degrees 
of freedom and thus constitute heat baths. As a conse- 
quence, phenomena such as energy relaxation, fluctua- 
tions, and dephasing influence the system dynamics. In 
the classical domain the corresponding theoretical de- 
scription is well developed and based on Langevin equa- 
tions or, equivalently, Fokker-Planck equations for phase 
space distributions [1| . The issue of dissipative quantum 
system has been initiated in the 1960s in the context of 
nuclear magnetic resonance [1, |1] and quantum optics [1] 
and since the 1980s has attracted substantial research 
in condensed matter physics 0, @. Most recently, it 
has regained considerable attention for quantum infor- 
mation processing, where noise appears as an undesir- 
able, but yet unavoidable effect [Tj]. In these systems the 
coupling between the relevant system and environmental 
degrees of freedom is designed to be as weak as possible. 
In contrast, finite dissipation can also be constructive 
and lead to qualitatively new processes such as stochas- 
tic resonance charge transfer in molecular structures 
: or ratchet induced transport in biomolecules [ll| ■ 
The regime of strong friction is classically known as the 
Smoluchowski limit .12.] and corresponds to a separa- 
tion of time scales between equilibration of momentum, 
which is fast, and equilibration of position, which is slow. 
The generalization to the low temperature quantum do- 
main, only a few years ago given in jlsj . shows that 
quantum fluctuations may appear at relatively elevated 
temperatures and substantially influence the dynamics. 
Since then this issue has triggered various applications 
for strongly condensed phase systems, see e.g. [l4| - [l9| . 

Theoretically, the description of quantum Brownian 
motion possesses an exact solution within the path in- 
tegral representation for the reduced density matrix 
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[S S HSl • This expression reveals that due to the non- 
Markovian nature of quantum mechanical fluctuations a 
simple equation of motion for the reduced density does 
in general not exist. Progress can be made in the weak 
coupling regime addressed above, where powerful master 
equations have been derived and successfully applied e.g. 
in quantum optical systems [1]. With a typical system 
frequency denoted by wq, the condition for this reduc- 
tion reads ^ huio,k^T, where 7 is a typical cou- 
pling strength between system and bath and T is the 
temperature. In the opposite range hq 3> hiOQ^k^T fric- 
tion dominates such that, roughly speaking, the typical 
linewidth induced by the environment exceeds the bare 
line separation as well as the thermal energy. In this deep 
quantum domain, named quantum Smoluchowski range 
(QSR), the reduced dynamics is nearly classical, but with 
a substantial impact of quantum fluctuations. Since fric- 
tion dominates, any approximate treatment must start 
from a formulation where the system-bath interaction is 
described non-pcrturbatively as e.g. in the path integral 
representation. It turns out that within this formulation 
a semiclassical type of approximation applies and the re- 
duced dynamics can equivalently be cast into an equation 
of motion for the marginal distribution in position, the 
so-called quantum Smoluchowski equation (QSE) |2l[. In 
leading order the quantum fluctuations in the QSE have 
been derived in [13| . 

What has not been done yet, is a systematic analy- 
sis of higher order corrections to the original QSE. In 
principle, this is a formidable task as it requires a sys- 
tematic expansion of the full path integral expression for 
the real-time dynamics. However, as already pointed out 
in , if a time evolution equation for the position dis- 
tribution P exists at all, it must have the form of a con- 
tinuity equation, i.e. P = CP. For the distribution in 
thermal equilibrium Pp one then has CPp = which 
can also be seen as an equation for C provided Pg is 
known. Higher order quantum corrections in C can thus 
be determined from a systematic approximation to the 
exact path integral expression in imaginary time of the 
reduced equilibrium density matrix. Additional dynam- 
ical corrections in C need then be analyzed in the time 
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window 1/7 <IC t <C 7/1^0 only. The corresponding ex- 
tended QSE covers the dynamics of the position distribu- 
tion for strong friction from high temperatures jhp <C 1 
to low temperatures jh/S ^ 1. 

An alternative approach, recently proposed for the 
high teniperature range 7^-/3 ^ 1 by Coffey and co- 
workers |22h24| , follows a similar strategy but is based on 
the thermal Wigner distribution for the uncoupled sys- 
tem. While this allows to obtain the universal leading 
quantum correction, we show here also, that it is not a 
consistent procedure to treat higher order quantum cor- 
rections. The latter ones carry information about the 
system-bath coupling, which is absent in the bare ther- 
mal distribution. 

The paper is organized as follows: In Sec. HI] a brief re- 
view of the path integral representation is given and the 
relevant regimes in parameter space are discussed qual- 
itatively. Then, in Sec. IIIII the generic case of systems 
with harmonic potential is studied for which the real- 
time dynamics can be solved exactly. This allows for 
a detailed analysis in the strong friction range at low 
temperatures. The extension to anharmonic systems is 
presented in Sec. IIVI where the perturbation theory in 
the inverse friction strength is applied to derive higher 
order quantum corrections to the Moyal coefhcients in 
the QSE. Interestingly, this type of semiclassical analy- 
sis also provides a formal solution to the problem which 
reveals the role of minimal action paths and deviations 
around them. As an application the escape rate out of a 
metastable well is calculated and shown to reproduce in 
higher order perturbation theory the exact result for the 
quantum enhancement factor above the crossover tem- 
perature. The paper concludes with a discussion of al- 
ternative approaches in Sec. |V] and a summary of the 
results in Sec. |Vll 



II. PRELIMINARIES 

A. Dissipative quantum dynamics 

Dissipative quantum systems are described with sys- 
tem-|-reservoir models [2^ , where the position g of a 
system with potential V{q) is bilinearly coupled to the 
positions Xa of environmental oscillators. Therefore, the 
Hamiltonian reads 
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(1) 



with the g^-dependent term in Hj added in order to 
avoid coupling-induced potential renormalizations. The 



dynamics of the full system described by a density oper- 
ator W{t) is then given by 



W{t) 



W^(0)e+^ 



Ht 



(2) 



The initial state W{Q) is obtained from the equilibrium 
density operator Wjs through the application of projec- 
tion operators acting on the Hilbert space of the sys- 
tem only. In contrast to the Feynman- Vernon theory, 
it therefore bears initial correlations between the sys- 
tem and the bath. We are, however, only interested in 
the reduced system described by the density operator 
p{t) = trB{VF(t)} and particularly focus on its position 
representation {q\p\q') = p{q,q'). Further, the analysis 
will be restricted to a class of initial preparations given 

by 



P{qi,q'i,t = 0) = P0{qi,q'i)Mqi,q'i) 



(3) 



with a two- variable preparation function X{qi^q[). Stem- 
ming from projection operators, it describes deviations 
from the reduced thermal equilibrium density pp = 
trB{exp(— /3i?)}/(ZZB) with the bath partition function 
Zb and a proper normalization for the system Z . The 
path integral approach allows for an exact elimination of 
the bath degrees of freedom in the position representa- 
tion. In doing so, one obtains 



piqf,q'f:t) 



dqtdq-J{qf,q'f,t,q^,q-)X{q„ql) . (4) 



The propagation function J{qf,q'^,t,qi,q'j) contains a 
threefold path integral over the system coordinates, two 
for the real-time propagation of the initial density ma- 
trix, one thermal path in imaginary time for the initial 
state. For further details we refer to the literature f5l.[20|. 
In time t, the real time paths run from qi and q[ to qf 
and q'j^ respectively. On the imaginary time axis, the 
thermal path runs from qi at time to q[ at time —ihf3, 
where /3 = {kBT)~^ . Besides the bare propagation, the 
integrand in the propagation function contains an influ- 
ence functional that keeps track of the interaction with 
the bath. The latter one is nonlocal in time and contains 
the damping kernel 



KiO) 



duo cosh [w(?i/3/2 — i6 
~ sinh(ft/?w/2) 



(5) 



which is proportional to the force-force correlation func- 
tion of the bath. It is determined by the spectral density 



2maUJc 



-S {U! — UJa) . 



(6) 



For purely real time arguments, this kernel reduces to 

K{t) = / — /(w) [coth{h(3u}/2) cos{ujt) - isin(wt)] , 
Jo 

(7) 



3 



and for imaginary times one finds K^—ir) = 
fi : 6{t) : — fc(r) with a time local contribution fi — 
{2/tt) J duj I (uj) / u! . The time nonlocal one has a represen- 
tation in terms of Matsubara frequencies Vn — 2T:n/h/3, 
i.e., 



fc(r) 



m 

hp 



(8) 



where 7(2) denotes the Laplace transform of the classical 
friction kernel 7(i). 

In the sequel, we consider Ohmic damping, i.e. a spec- 
tral density of the form /(w) = 77170; and employ a Drude 
regularization, I{uj) = mjojui'^/ (w^ + w^) with a high fre- 
quency cut-off frequency Wc, whenever needed. 



B. Regimes in parameter space 

While the path integral formulation provides with (|4]) 
a formally exact expression for the reduced density ma- 
trix, its explicit evaluation is very demanding. Analytical 
results have been derived for harmonic systems, but in 
general even numerical approaches are limited to specific 
models or regimes in parameter space. It has further 
been shown that an equivalent "simple", i.e. tractable, 
equation of motion for the reduced density does not ex- 
ist (see e.g. [1^). The reason for these complications is 
the non-Markovian nature of quantum Brownian motion 
as determined by the time nonlocal damping kernel ([7|. 
In this situation progress can be made at least in a per- 
turbative sense in certain ranges of parameter space. In 
this section we discuss the two complementary domains 
of weak and strong friction, respectively. 

For this purpose, a typical damping strength in the 
long time limit is introduced as 



7 = 7(0) = lim 



lieu) 



w->0 muj 



(9) 



such that both for Ohmic and for Drude damping we have 
7 = 7. A second relevant frequency scale of the bath 
is temperature, i.e. Vi = 2'K/h/3, while for the isolated 
system we assume a typical frequency ujq. In fact, these 
three scales define qualitatively the nature of the reduced 
dynamics. 

In case of weak friction j/ujo <C 1, the typical relax- 
ation time towards thermal equilibrium is of the order of 
tr oc 1/7 and thus much larger than any time scale of 
the bare system dynamics. As far as we are interested in 
phenomena like decoherence and dephasing, the damp- 
ing kernel ^ becomes time-local on a coarse grained 
time scale t ^ hfH provided that jh/S <C 1. This latter 
condition is the basic assumption for all types of master 
equations which have been used in fields like quantum 
optics, nuclear magnetic resonance, and most recently in 
quantum information processing. 

In the opposite domain of strong friction ^/ujq 3> 1 
again a time scale separation exists, since the relaxation 



of position happens to occur for times on the order of 
tr OC j/ujg, while the momentum relaxes on the time 
scale 1/7. Accordingly, a coarse graining of time leading 
to a time-local kernel ([7]) is possible if 



11 7 
^^'T' — < 



7 



(10) 



On the one hand this latter relation comprises the clas- 
sical regime <^ 1 corresponding for strong friction 
to LjQhp <^ 1, and on the other hand the deep quantum 
regime 7/1/? » 1. 

In the classical domain of strong friction the Fokker- 
Planck equation for the distribution in full phase-space 
reduces to the famous Smoluchowski equation for the 
marginal distribution in position P{q,t) = {q\pit)\q) 
[IlUJj namely. 



P{q,t)^—d, 
7717 



V'{q) + ^d. 



P{q,t). (11) 



The quantum regime (Quantum Smoluchowski regime, 
QSR) has gained much attention only recently after a 
quantum Smoluchowski equation (QSE) has been derived 
from the path integral expression in 13]. The aim of this 
paper is to systematically derive higher order corrections 
to the original QSE in order to better understand its 
limitations, but also to obtain improved results. 



III. HARMONIC OSCILLATOR 

As mentioned above, harmonic systems V{q) = 
allow for an exact solution of the path inte- 
gral expression (jj]). They may thus serve as models to 
analyze the QSR and the existence of a QSE in detail. 



A. Time-dependent density matrix and time 
evolution equation 

According to [1^ the propagating function for the re- 
duced dynamics can be calculated explicitly and the re- 
duced density matrix ^ follows from 



p{rf,Xf,t) 



47r|A(i)|y2^ 



dxidnX{ri,Xi 



X e 



(12) 



where we introduced sum and difference coordinates r = 
{q + q')/2 and x = q — q', respectively. Here, the minimal 
action S(-) is of Gaussian form in the coordinates. It can 
be expressed in terms of the symmetric and antisymmet- 
ric part of the position autocorrelation function 

S{t)=Re{{q{t)q)p}, A{t)^lm{{q{t)q)p} (13) 

which are determined by the natural frequencies 7/2 ± 



Further, the mean square of position in 
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thermal equilibrium is found to read 



1 ^ 1 



mP ^ 1/2 + 1 1/„ |7 + 

n— — oo 'til' u 



(14) 



In general, progress can now only be made upon spec- 
ifying the preparation function in (jl2p explicitly. How- 
ever, simplifications arise in the strong friction domain 
considered here. Namely, for large friction 7/a;o ^ 1 and 
on the coarse grained time scale [cf. (fTOj) ] 



-, — < t 

7 '^c 



only the frequency 



(15) 



(16) 



is relevant. In this parameter regime, the position auto- 
correlation functions reduce to 



4m(7/2-r2) 

m = - 



cot {hl3n/2)e 



-Qt 



-Qt 



(17) 



so that e.g. A is of order I/7 or smaller since 7/2 — £7 « 
7/2. Further, in leading order the variance in position 
reduces to its classical value (q^) w l/m/J^Q. Then, a 
straightforward analysis reveals that the Gaussi an fa c- 
tor in E restricts Xi to be at most of order A/ y/ (g^) ~ 

l/vw(?) ■ Accordingly, the x^-dependence of the prepa- 
ration function is only probed on this short length scale 
and we may set X{xi,ri) « A(0,ri) in ([T^ provided the 
initial preparation is sufficiently smooth in Xi. Equiva- 
lently, the initial preparation must be restricted to mo- 
menta of order 7-^/ (g^) or smaller. This, however, is in 
complete accordance with the Smoluchowski limit which 
requires equilibration of momentum on sufficiently short 
time resp. length scales. 

After performing now the -integration and keeping 
the dominant terms only one arrives at 



p{rf,Xf,t) 



\/27rA 



X exp ■ 





\ S{t)] 




1 2A 







where 



A 



1 



sHt) 



(18) 



(19) 



and Ppiri) — pplji^Xi = 0) denotes the diagonal part of 
the reduced thermal distribution. The equilibrium vari- 
ance in momentum is given by 



+ 00 

E 

n— — oo 



(20) 



This sum must be regularized by introducing a cut- 
off frequency Wc as discussed above. In the inter- 
esting range 7?i/3 ^ 1 we then easily find (p^) ?a 
(m?i7/7r)log(cJc/7), while in the complementary high 
temperature regime 7?i/3 <§; 1 one regains the classical 
results (p2) « (TO//3)[l + (n/37/7r)ln(n/3wc/27r)]. The lead- 
ing corrections to the result (IT51) consist of a phase factor 
exp[(i/?i)x/(5p(r/, r^)] with 



5p(r/, 



A(t) m{q'^ 
A(t)~A' 



(q2 



(21) 



For strong friction one has 51 sa '^q/i so that A/ A « 
l^q/j- The order of magnitude of the factor in square 
brackets is given by the width of the Gaussian in ([T8|) as 
VA and that of by the Gaussian prefactor as 1/ ^/(jp). 
Hence, the above phase factor produces corrections in the 
exponential on the order of 1/^/7^ (p^) or smaller. These 
corrections depend on the cut-off parameter Wc and can 
be made arbitrarily small. Other corrections to ^TE\\ are 
at most of order exjp{—jt) and are negligible on the coarse 
grained time scale. 

Now, a Wigner transform of (fT8)l immediately 
gives rise to a momentum distribution proportional to 
exp(— p2/2(p2)) as expected. For the relevant diagonal 
part P(q,t) = p{rf,Xf = 0,t) the dynamics can then be 
cast into an equation of motion for its time evolution 28 1 



Hq.t) 



dg[Di{q) + dgD2]P{q,t) (22) 



,2^9 



with a drift coefficient Di{q) = muj'^q and a posi- 
tion independent diffusion term D2 = raLUQ{q^). To 
leading order one regains from the above expression in 
the QSR the classical Smoluchowski equation (|lip with 
_D2,ci = 1//? and ^/muj"^ w Xjm^. A systematic ex- 
pansion around this result will be given in Sec. IIII Dl 
Moreover, for systems driven externally by time depen- 
dent forces f{q,t) the analysis goes through accordingly 
by replacing mujQq — >■ mujQq-\- f{q, t) provided the typical 
time scale for the driving sufficiently exceeds the scales 
1/7, hf3. This result has also been derived from a Quan- 
tum Langevin equation in [29| . 



B. Initial correlations 

The pioneering work on quantum Brownian motion 
by Feynman and Vernon uses a factorized initial state, 
where the initial density of the whole system takes the 
form W{Q) = ps{0) exp{—f3HB) /Zb ■ Initial correlations 
between system and reservoir are thus absent. Asymp- 
totically, for very long times the reduced density matrix 
approaches the true reduced thermal equilibrium distri- 
bution as well, but for intermediate times factorizing and 
non-factorizing initial conditions lead in general to differ- 
ent stochastic processes. It is known that for high tem- 
peratures/weak friction both processes coincide after a 
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transient period of time which is short compared to the 
relevant dynamics [1, [1]. Question here is if this holds 
true also for low temperatures/strong friction so that the 
QSE can also be derived within the somewhat sim- 
pler factorizing formulation. 

Following the lines described in the previous section 
one gains for factorized initial states a result which looks 
similar to (|T8)) with the replacements X{0,ri)Pp{ri) — >■ 
ps{n,t = 0), S{t)/{q^) -2mA{t)/h, and A ^ A = 
AA'^{p'^)/h'^ + A. In contrast to p8)) . however, an addi- 
tional phase factor exp[{i/ h)x fr fUiA/ A] appears, which 
cannot be neglected. Namely, while the order of mag- 
nitude of the corresponding phase factor (j2ip for non- 
factorized initial states is completely determined by the 
widths of the Gaussians in Xf and r/ — riS/ (g^), respec- 
tively, this is not the case here. The order of magnitude 
of this term depends via r/ also on the distribution of 
initial values for r^. Hence, the density matrix does not 
reduce to a simple product of thermal momentum distri- 
bution and time dependent position distribution. 

If one ignores this complication and focuses on the di- 
agonal part of the density p{q, Xf = 0) only, a time evo- 
lution equation of the form (j22p is found, however, with 
a diffusion coefficient which directly depends on the mo- 
mentum variance via {p^) A{tY / {q^) . While this term is 
negligible in the classical domain (it is then of order 1 /7^) 
and for times t ^ 7/1^0 ' provides an essential contri- 
bution for 7?i/3 ^ 1 on the coarse grained time scale ([T5|) . 
The result is a Smoluchowski type of equation which con- 
tains via a time dependent drift coefficient the momen- 
tum variance. This analysis reveals that the formulation 
based on factorized initial states does not lead to an ac- 
ceptable Smoluchowski description in the deep quantum 
domain. The reason for this behavior is that for strong 
friction the time scale on which initial correlations are 
established is identical to the relaxation time of the full 
system, i.e. equilibration of position. 



C. Equilibrium distribution and current operator 

Since the treatment of the dynamics laid out in 
Sec. nil Al requires ergodicity, a time evolution equation 
for the reduced density must describe the relaxation of 
the system to thermal equilibrium. As we will see here, 
this allows for an alternative derivation of quantum cor- 
rections to the Smoluchowski equation starting from the 
equilibrium density matrix of the open system. 

The time evolution equation (|22l) has the structure of 
a continuity equation for the marginal probability distri- 
bution in position. Accordingly, the thermal distribution 
Pp of the harmonic system is determined by a vanishing 
current J = LPf^ici) — with a current operator being 
defined as 

C^DM+dqD2 (23) 
with Moyal coefficients Di and Z?2- Note that only for 



systems with at most quadratic potentials D2 is a con- 
stant. 

The idea is now to consider (I23p as an ansatz with yet 
unknown Moyal coefficients. The thermal distribution 
has thus to be of the form 



1 p-tpiq) 



D2 



(24) 



with Zq being the partition function of the oscillator. 

On the other hand, the reduced thermal density matrix 
for a general system can be represented as a path integral 
in imaginary time, namely. 



1 



V[q] e- 



f[9l/?i 



(25) 



where q(r) connects in the time interval h/3 end-points 
q{0) = q' with q{hf3) — q. The effective action contains 
the bare euclidian action of the system and the influence 
functional in imaginary time determined by the kernel 
([5]). Z denotes a proper normalization, which in case of 
a harmonic system coincides with the partition function. 
For quadratic potentials this expression can be evaluated 
exactly to read 



1 



exp 



2(g2 



2^2 



(26) 



with the variances in position and momentum, respec- 
tively, given in (ITil) and (PH)) . respectively. Upon com- 
paring Pp with the explicit result for Pp{q) = pp{q,0) 
one can determine Di and D2. As expected, for the har- 
monic oscillator one regains the results specified already 
in (|22l). 

This alternative procedure reveals an important fea- 
ture. We can distinguish two types of corrections to the 
classical Smoluchowski equation: Those in the current 
operator are completely determined by the equilibrium 
properties of the system and are influenced by quantum 
fluctuations; the other ones are dynamical in origin and 
appear merely as an overall factor in front oldqC [cf. (p2t ] 
completely determined by classical dynamics. Hence, as 
long as we are interested in the role of quantum fluc- 
tuations in the QSE for anharmonic systems, it may be 
justifled to perturbatively calculate the Moyal coefficients 
from the corresponding thermal distribution and assume 
that quantum corrections in the dynamics are relevant 
only in higher orders of perturbation theory. In fact, the 
analysis presented first in [l^ proves this scheme to be 
correct at least for the dominating quantum fiuctuations. 
To prove that it applies also to higher order contributions 
requires to consider the reduced dynamics only within the 
time window I/7, h(i <^ l/^vi- 



D. Perturbation theory for strong friction 



So far the QSE obtained in (1221) is exact on the coarse 
grained time scale. To lay the basis for a perturbative 
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treatment in case of anharmonic potentials, we study 
in the sequel a systematic expansion around this exact 
expression in terms of the small parameter wo/7- The 
corresponding results are thus valid for all temperatures 
provided damping is strong. For this purpose, one writes 



where the coefficients serve as expansion parameters. 
We emphasize that this approach covers the high tem- 
perature regime "fh/3 ^ 1 as well as the challenging deep 
quantum domain 7/1/3 3> 1. In contrast, an alternative 
formulation developed recently [13] is restricted to the 
former range. For further discussions we refer to Sec. [V] 



mp 



Tl = l 



(27) 



where ((Z^)ci = ^/iPmujQ) is the classical variance in po- 
sition. For strong friction 7/wo ^ 1 one then has 



A 



ml3 ^ ^ ^ 



n=l 



(28) 



x=0 



The coefficients are independent of system properties 
and can be expressed with the help of the polygamma 
functions. In lowest orders they read 



An 



Ai 



2 



^-(1 + - 



Ce 



2v 

7 



* 1 



Ce 



(29) 



with Ce being the Euler-Mascheroni constant and "^'^^^ 
the /cth derivative of the digamma function 5*. In 
particular, one has in the high temperature domain 
7^/3 <C 1 the approximation Ao ~ ?i^/3/12m, while in 
the quantum regime 7/1/3 ^ 1 the expression Aq ~ 
(/i/7rm7)ln(7/i/3/27r) reveals that quantum fluctuations 
are substantial and not algebraically small. The coeffi- 
cients A^,/i > 1 are on order of ln(7)/7^'^"'"^ and thus 
serve as semiclassical expansion parameters. We note 
that their dependence on Planck's constant is highly non- 
algebraic. This way, one obtains the leading quantum 
corrections in (|22]) as 



-d„ 



{l + mulp\o)dq P{q,t), 
(30) 

where the factor in the diffusion coefficient may in 
this approximation also be written in the form 1/(1 — 



A. Perturbation theory 

In the spirit of the previous section we consider the 
equilibrium density (j25p and obtain for its dominating 
contributions a result of the form 



A'r.x) = ^exp 



n 



lir, x) 



Fmair,x), (31) 



with the minimal action SE,ma calculated from the min- 
imal action path qrna{T)- It is given by 

mq,na - V'{q„^a) - / CLt' k{T - T')q„-,a{T') = . (32) 
"'0 

subject to the boundary conditions ^(0) ~ r + x/2 and 
q{fif5) — r — x/2. Fluctuations around the minimal action 
path are accounted for by 



^ma •^) 



y{m=o 



2?[y]e-5fe/^''^-J/M(£i')M (33) 



/y(0)=0 

containing the second variational order operator 



[-m92 + V"iqrna)] y{T) 

+ / dr' k{T - T')y{T') . 
Jo 



(34) 



In case that there are several minimal action orbits one 
has to sum over them in This general expression 

must now be specified to the situation in the QSR. The 
underlying picture there is this: The strong friction sup- 
presses off-diagonal elements of the density matrix mean- 
ing that the minimal action path does not deviate much 
from the mean r of its initial and final position. Thus, 
it suffices to restrict the analysis to the diagonal part 

Ppiq) = Pi3{q,o), i.e.. 



Pp{q) = ^F^a{q,0)e-^^''^/^ 



(35) 



with S{q) = Sma{q,0)- The minimal action path and 
deviations from it are then treated perturbatively in the 
spirit of a self-consistent approximation. For this purpose 
it is convenient to switch to Fourier space 



IV. ANHARMONIC SYSTEMS 

Departing from the harmonic oscillator potential, we 
have to resort to approximation schemes in order to solve 
the path integral expressions. In the strong friction limit 
a type of semiclassical or saddle point treatment applies 



leading us from ([5^ to 

^lln +7Wn\qn + Vn/m ■ 



b. 



(36) 



(37) 
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Here, 



Jo 



(38) 



and the term on the right hand side accounts for the fact 
that the Fourier series periodically continues the path be- 
yond the interval (0, hfS) producing singularities at times 
T = nhp, n = 0, ±1, ±2, . . . with b = qma{fiP) - gma(O). 

The strong friction approximation requires systems 
with sufRciently smooth potentials so that locally the ex- 
pansion holds 



V'[qrna{T)] (q) + V" (q) [qmair) - q] 
V"'(a) 9 



(39) 



From (|36| and ((37)) one observes that the qn,n ^ are 
suppressed at least by factors of order I/7, while q^ must 
be of order 1 due to 9ma(0) — qma{fi(3) — q- Conse- 
quently, deviations Sq — qma(j)—q are suppressed by fric- 
tion as well. Introducing an anharmonicity length scale 
I the above expansion thus serves as a starting point for 
a systematic perturbation theory provided that \5q\ <C I. 

With the above strategy we obtain an explicit expres- 
sion for the thermal distribution in the overdamped limit. 
To determine with this result quantum contributions in 
the QSE, we follow Sec. UTTCland look for a QSE in the 
form P — dq^C P with a dynamical factor ^ and a cur- 
rent operator C ^ Di{q) + dqD2{q)- Now, if one writes 
P/i in the form 



1 Q—^ii) 



(40) 



a vanishing equilibrium current [Di{q) + dD2{q)] Pjs = 
is guaranteed in all orders of perturbation theory. Upon 
comparing the semiclassical expression ([55]) with (UHl) the 
most obvious choice is to identify ipiq) — S{q)/h and 
D2 oc l/Fmaiq^O) with a proportionality constant fixed 
by the high temperature limit. In this range S{q) — >■ 
Hf3V{q), while the fluctuation factor tends towards the 
free particle result Fmaiq, 0) ■\/m/27rfi?/3. Hence, with 
the scaled fluctuation factor 



F{q) ^ F^a{q,0) 



(41) 



we have the formal expressions for the Moyal coefficients 



pF(q) 



PF(q) h dq 



which classically reproduce Di{q) — > V'{q) and D2 — 

The quantum mechanical current operator is not 
uniquely fixed by the above choice for the Moyal coef- 
ficients though. In fact, any operator of the form 



obeys CPp = and written as C 
define new Moyal coefficients 



Di + dqD2 one can 



b^{q) - d{q)D^{q) - d'{q)D2{q) , ^2(9) = d{q)D2{q) . 

(44) 

Thus, it seems that there is a whole class of Moyal co- 
efficients, where each pair follows from the representa- 
tive one (|42|) by a proper function d{q) obeying d — > 1 
in the classical limit. By construction all these current 
operators produce a vanishing current for the thermal 
equilibrium ([35]) . However, based on the previous sec- 
tion the current operator for a harmonic system is known 
exactly. This information can be used as an additional 
constraint to fix the correct d{q) as we will show in the 
sequel. Eventually, dynamical corrections ^ are obtained 
via an evaluation of the reduced dynamics Q within the 
time window l/j,hf3 t ^/ujq. 

One remark is in order here: While the representa- 
tion of the current operator as Di + dqD2 with the co- 
efficients pi)) may be advantageous for physical inter- 
pretations, it has one major drawback compared to P5)) 
though. Namely, since ip = J dyDi/D2 = S/h— \n{d) a 
numerical calculation of ip based solely on the knowledge 
of the coefficients D1/2 gives an exponent for the equilib- 
rium distribution which carries in addition to the correct 
action term lower order terms contained in d. As d is 
not known explicitly in this representation the numeri- 
cally obtained ip may give rise to spurious equilibrium 
currents (see e.g. |14l|). Hence, particularly in numeri- 
cal applications the representation (|43l) is superior. We 
will illustrate this by discussing explicit expressions in 
Sec. UVCl 



B. Leading order 

In leading order we put in ([39]) V [q™, 
with Vn w V {q)h(3Sr,,o we find = 
n 7^ 0, one has 



,(0) 



bio) 



a] — V'{q) so that 
hpV'{q)/m. For 



(45) 



C = d{q) [D^{q) + dqD2{q)] 



(43) 



and go is determined via p6l) from (7^0 (0) = 9- This leads 
to q^^ /hp = q — fo^^-* rnXoh. The deviation Sq can there- 
fore be estimated to be on the order of {h(3/j)\n{'-fh/3) for 
jhf3 ^ 1 and on the order of {hf3)'^ for <C 1, thus 
justifying the perturbative treatment. 
The corresponding minimal action reads 

S(o\q)^hmq)-^V'iqr. (46) 

The fluctuation integral ([55]) is calculated by replacing 
V^"[9ma] — V"{q). The boundary conditions are most 
conveniently taken into account by introducing (5[y(0)] = 
S[J2n Vn] in the path integral and representing the 5- 
function as an integral over an auxiliary variable. This 
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way, one gets f ("^(g) = [1 - l3V"{q)\o\/ 13 and the ther- 
mal distribution in the form pj^"^ ~ F^^^ Pg'^'''/^ with the 
quantum correction 



= [l-/3r'(g)Ao]e(^'/2)y'(,)^Ao 



(47) 



and the unnormalized classical distribution P, 



id) 



exp(-/3y) (see also fig. O in Sec. UVE)) . Now, fi-om (02) 
we immediately derive 



1 



1 



p l-/31/"(g)Ao 



(48) 



The dynamical factor follows from the path integral cal- 
culation as in and reads as in the classical case 
^ = l/m7. Note that in this order of perturbation the- 
ory Df\q) w [1 + l3V"{q}\Q[J3 so that we regain the 



QSE already obtained in [131, l3(l|. For a purely harmonic 
system the result reduces to (|30)) . 



C. Local harmonic approximation 

In next order the full potential is approximated around 
the endpoints g as a harmonic oscillator meaning that 
in ([M)) all terms beyond the V"{q) contribution are ne- 
glected. From ([57]) we obtain for n ^ 



,(1) 



fed 



(49) 



In analogy to the preceeding section, we define a position 
dependent function 



1 



(50) 



which can be expanded according to p8| with Wq re- 
placed by V"{q)/m. Summation over all Fourier com- 
ponents of the minimal action path leads together with 
the boundary condition to q'^l^ = h/Sq — pb^^^mA{q). The 
velocity jump at the endpoint follows as 



hp 



V'(q) 



l + PV"{q)K{qy 



(51) 



Expanding the minimal effective action up to first order 
in 5q = qrnair) - q yields 



S('Hq)^'^qb^'^+hp [v{q)-^V'{q) 



V'{q) qV'iq) 



hp 



dr [qma{T) - q] 



--hl3V{q) - h0- 



,V'^{q) 



A(<Z) 



l + l3V"{q)K{q) 



(52) 



The fluctuation path integral F{q) is treated accordingly 
with the second derivative of the potential expanded in 



the same way as discussed above. In the local oscillator 
approximation, first and higher order terms in q^a (''") ~ 9 
are truncated consistently. Defining u„ = [i^^ + 7|i/,i| + 
y'(<7)/m]~i, one has for the scaled factor 



+ 00 



+ 00 



-1/2 



Uo 



exp {-(m/3/2) E^>i A^_i [V"{q)/mr / ^,\] 



^l + PV"iq)Aiq) 



(53) 



Now, the expressions ((52)) and ((53)) combine to provide 
the position distribution in thermal equilibrium for a 
strongly overdamped system ((35)) . Since this result is 
interesting in its own as it applies to all systems with 
sufficiently smooth anharmonic potentials, it is worth to 
look for quantum corrections in more detail. Upon ex- 
panding the function A{q) similarly as in ((^5)). quantum 
mechanical contributions to the classical distribution can 
be systematically derived. Taking into account terms 
up to order in Aq wc have for the quantum correction 



y(l) ^^pW/p^cO 



X exp 



1 - pV'iq) 
P'V'iq) 



3pV"iq) 



2 r 



Ao + 



4m 
V"{q) 



\^-PV"{q)\l 



(54) 



In this expression Ai and Ag dependent terms have been 
kept so that it applies to the entire strong friction range 
with 7/wg S> hp. In the high temperature regime 
7/1/3 ^ 1 one has Ai ^ mp\^, while for low temper- 
atures 7/1/3 ^ 1 the inverse relation holds and Y^^'^ be- 
comes a pure perturbation series in powers of Aq. In fig. [3] 
in Sec. IIVE) the impact of these quantum fluctuations is 
exemplified. 

Before we proceed, we mention that the approach dis- 
cussed so far can also be extended to include off-diagonal 
elements of the density matrix. The straightforward cal- 
culation leads to 



"(^'(r,a;)=p(^'(r)e-"('-)-'/'' 



where 



' n=l 



V"{q) 



(55) 



(56) 



corresponds to the local momentum variance. This func- 
tion can be expanded according to 



^=1 ^■ 



(57) 
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For purely ohmic damping 770 diverges and must be reg- 
ularized by e.g. a Drude cut-off 0;^ ^ 7 yielding 



mh'j 



V 2v 

27 



(58) 

In the ohmic case the first coefficient stays finite and 
reads 



(59) 



Now, that the thermal distribution is at hand, the 
Moyal coefficients follow according to (H^ with F^^^ 
specified in (l53l) and with the action as in (|52|) . To 
obtain the function d{q) in the ansatz ()43p . one ob- 
serves that for a harmonic system the exponential in 
F^^) is up to a temperature dependent factor identical 
with its partition function Zq, while the nominator is 
y/W)lWhi- Thus, the known result ^ is regained if 
d^^-' = F^^^ (1 + Pmbj'f^K). The generalization to anhar- 
monic potentials leads then to 



d(l)(g)=F(l)(q) [l+pV"{q)A{q)] 



(60) 



Of course, in the classical high temperature domain one 
has d!^^^ — 1. As discussed around (H^ . the current op- 
erator C = d[Di + dqD2\ with the known d can also be 
cast into the standard form Di + dqD2- The result is 



dq 



D 



(1) 



- [l-f/3r'A(g)] 



(61) 



With these expressions the current operator in the QSE 
is completely determined. For practical applications it 
is more convenient though to systematically expand the 
above results to the desired order in the inverse friction 
strength. For this purpose we recall that formally the 
^inl^ > 1 are of order ln(7)/7^^+^, while Aq is of order 
ln(7)/7. In particular, in the QSR range ^h[i 2> 1 one 
obtains 



D['>iq) « V'{q) \l + ^Xl [V"{qr + 2V' {q)V"' (q)] 



D^\q)^i 



P 1 - PV"{q)\o + lP^V"{qf\l 



S'Hq)^l~^^V"{qf\l. 



(62) 



Quantum fluctuations appear in the effective drift coeffi- 
cient 

D^'\<l) « V'{q) + ^\lV"\q) [V"{q)+f3V\qf] (63) 

via the anharmonicity of the potential, while in d''2^ sa 
[1 -I- (3V"\q\/(5 only harmonic properties are contained. 
Corrections here are of order Ai ^ Aq. 



At this point we come back to our discussion at the 
end of Sec. IIV AI We consider e.g. a periodic potential 
for which one immediately sees that = S/h is periodic 
as well in each order of perturbation theory for the ther- 
mal distribution. For ip = J D1/D2 which must be inte- 
grated numerically this is by no means obvious. To con- 
sistently neglect in the numerical ip contributions higher 
than of order Aq is impossible which may thus lead to fi- 
nite currents in thermal equilibrium. This illustrates why 
it is advantageous particularly in numerical calculations 
to work with the coefficients ([62| in the current operator. 

In principle, the thermodynamic analysis presented 
above must now be supported by a dynamical calcula- 
tion of the reduced density (HJ to reveal whether quan- 
tum corrections in the dynamical factor ^ are rele- 
vant. This can indeed be done within the time window 
h/3,l/j ^ t ^ 'y/i^o- Effectively, the calculation goes 
through along the lines described for harmonic systems 
with ujq —i' V"{q)/m. As a result dynamical deviations 
appear in a form similar as in (|22p and are thus com- 
pletely classical; one has £, = 1/(7717) + 0{V" /y^). The 
QSE in this order of perturbation theory is thus obtained 
as 



P{q,t) = —dqC<^'^ P{q,t) 



■m-f 



(64) 



with the current operator defined either in the form (|43p 
with the coefficients (H^ and (|60p or in the standard form 
with (1611). 



D. Application: Quantum escape rate 

A non-trivial case to prove the consistency of the 
above QSE is also to calculate the escape rate out of a 
metastable well and to compare with the semiclassically 
exact result |3l|. The latter one can be derived e.g. 
within the ImF approach. Above the crossover temper- 
ature one finds [s^l 



2tt ujii 



(65) 



where ojq = V"{0)/m and ujf — \V"{qb)/m\ denote fre- 
quencies for small oscillations around the minimum of 
the well at q = and around the barrier top at g = g^, 
respectively. For the quantum enhancement factor one 
has 



+°° 2 I 2 
_ -r-r v^+UJ^ + ^nl 



(66) 



and ijjR = — 7/2 -f yjiol + 7^/4 is the classical Grote- 
Hynes frequency. An alternative representation for fq 
which is more convenient to compare with strong friction 
results is given by 



fq = exp 



9 2^ n\ 



(^o" + (-l) 



At-l-l 



(67) 
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0.02 



w 0.01 h 



0.0 =■ 




-2-10123 
ln(7/u;o) 

FIG. 1: (Color online) Escape rates according to the QSE 
compared to the semiclassical exact rate vs. friction strength 
for a metastable potential with ij;o = in the high temper- 
ature range. The solid line depicts the leading order result 
[Ao-term in (|67p ] for temperatures ijjohfi = 0.2, 0.4, 0.6; dif- 
ferences between the three data sets are indistinguishable on 
this scale. Also shown are results obtained according to [2^ 
with a friction independent quantum enhancement factor for 
uiohp = 0.2 (dotted), 0.4 (short-dashed), 0.6 (long-dashed); 
see text for details. 
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FIG. 2; (Color online) Same as in fig. [T]but in the low tem- 
perature range uohp — 5 and only for the QSE. The leading 
order approximation (dashed line) and the next order approx- 
imation including terms of order A2 (solid) are shown. 



rameters leads to 



7mJ 
DM 



-"4'iq} 



dy 



Qi>{v) 



(68) 



Thus, with the normalization Z = Z^^eii of the harmonic 
oscillator in the well, we gain the escape rate in the form 



well 



7m 



9b 



dq 



dy 



D2{q) 



d{y) 



(69) 
(70) 



Upon evaluating the integrals in saddle point approxima- 
tion which is justified for sufficiently high barriers, one 
finds 



QSE 



27r7 



(71) 



This is a remarkable result as it reveals that the QSE re- 
produces the quantum enhancement factor (j66p exactly. 
The Grote-Hynes frequency corresponds to classical dy- 
namical corrections and for strong friction reduces to its 
known leading order approximation ljji « I/7. Hence, 
our thermodynamic procedure to derive quantum correc- 
tions in the QSE has also been proven a posteriori. In the 
following subsection we will show that we can extend this 
method even beyond the local harmonic approximation. 

The above findings are illustrated in figs. [H [51 where 
the ratio Vqse/^ is depicted for different orders of per- 
turbation theory. In the high temperature range (fig. [1]) 
already the leading order correction of the QSE gives ba- 
sically the exact result for all damping strength. In con- 
trast, in the approach recently proposed in |23{ the rate 
carries the enhancement factor for 7 = 0. This leads 
to increasing deviations for increasing damping strength 
and/or decreasing temperature and thus reveals the rel- 
evance of friction even at elevated temperatures (for fur- 
ther discussions we refer to Sec. |V|. For low temper- 
atures, see fig. [21 the QSE gives the correct result for 
the quantum fluctuations already for moderate friction 
strengths. Deviations remain small meaning that the 
perturbative expansion (|57)) quickly converges. Note that 
for potentials with = loj, all odd order contributions 
cancel. 



E. Higher order corrections 



The above expression relies on a local harmonic ap- 
proximation and the anharmonicity of the potential mat- 
ters only in so far as it leads to a variable curvature 
when one moves from the barrier top towards the well 
minimum. This situation thus perfectly fits to the type 
of perturbative treatment outlined in the previous sec- 
tion with a position dependent diffusion term. Station- 
ary non-equilibrium solutions of the QSE (I64p follow from 
C^-^'Tst = — J with a constant flux J. Variation of pa- 



We now take into account terms up to second order in 
the expansion pop . i.e. in the deviations 6q = qma{T) — q 
in the equation of motion for the minimal action path 
(|37p . Accordingly, one has for the zeroth Fourier compo- 
nent 



«f =h/3V'{q) + V"{q) iqo-h/3q) 



V'"{q) 



2E 



hp 



(72) 
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while for n 7^ coupling terms between different modes 
arise 



V"{q)qn + V"'{q)[^^-q]q,, 



E 



qa_ 

hp 

V'"{q) 



hp 



to first order in 5q leads to an integrand of the form 
exp [— (cr^^^ + ^tr)] , where the local harmonic part in 
the exponent ct'^^^ is given by replacing qma{T) q 
in V"{qma) as in the previous section. The next or- 
der term contains the third derivative of the potential 
(73) '5a = {V"'lh) J^^ dTy\T)[q{T) - q] and reads 



Since these coupling terms can be estimated to be of or- 
der V'^{q)V"' {q)\i, they are negligible against the other 
contributions for 7/1/? 3> 1. Hence, the leading impact 
of anharmonicities in the potential can still be treated 
analytically. Note that the term with \qm\'^ in the com- 
(2) 

ponent Uq is of the same order as the one with qmqn-m 
in Vn , but appears in (I37|) without friction and leads 
thus to a contribution which is larger roughly by a fac- 
tor 7 compared to the former one. In complete analogy 
to our treatment of the minimal action path in the local 
oscillator approximation, wc then find 



,(2) 



6(2 



'^n + lWn\+A{q)/m 



9 °° 1 
m n ^ — ^ 



"^^ ^1 + 7 \vn\ + A{q)/m 



, n^O (74) 



(75) 



Here, one observes the same structure as in the pre- 
vious section with the substitution V"{q) A{q) = 
V"{q) - F"'(g)m&(2)A(2)/S. Since A^^) is given by ^ 
only implicitly, we employ a linearization 



A(2) = A(g) + M , 6(2)^^(1)^^ 



(76) 



with A{q) as defined in (|50p to obtain in leading order 
SA w ^l3V"'iq)V\q)\oXi/m. Terms of order A{qf - Ag 
and smaller are neglected here, which also implies omis- 
sion of contributions of order A{q)6A ~ AgAi ^ Aq. In 
this approximation, the correction of the velocity jump 
reads 



mSb^- /3V"{q)mb'^^^SA 



PV"'{q) 

2h 



^6(l)A(g)' 



V"'{q)mbW^ 
2h 



Ai [l + f3V"{q)A{q)] 



(77) 



The action 5^2) — S'(i) -|- SS beyond the local harmonic 
approximation reads 



6S{q) = 



hP^V'iqfV"'{q) 
2m 



2^'('Z) - qV"{q) 



AqAi 



(78) 



which is of the same order of magnitude as 5 A. Formally, 
this correction is of order ln2(7)/(7'*^) with typical anhar- 
monicity length /. 

We now turn to the contribution of the fluctuations 
around the minimal action path Fma and discuss cor- 
responding anharmonic corrections. In the path inte- 
gral (j33p an expansion of the second variational term up 



5a- 



V"'{q) 

h 

, V"'{q) 



90 

hp-"^ 



E ■ 

m— — 00 



hp 



E 



{hpY 



-qn ■ 



(79) 



Here, again coupling terms between different Fourier 
modes appear and lead to a Gaussian path integral which 
is only tractable numerically. To make progress, we es- 
timate the various terms in (|79l) systematically. For this 
purpose, we recall that as far as orders of magnitudes are 
concerned we have {qo/hp~q) ^ Aq , qn ^ 7^ and 

Un ~ 1/^/7'"- 7^ 0, qQ,yQ ^ 1- The first sum in 6a pro- 
duces two types of contributions, namely, {qo/hp—q)y'^ ~ 
Ao and {qa/hp - q)yl^o ^ A0/7 ln(7)/72. The sec- 
ond sum gives rise to two types of terms as well, one 
with yo J2nyto ynq~n and one where the sum contains 
only ymy-m-nqn with TO,n, m -f rt 7^ 0. This latter 
part is estimated to be of order I/72 and thus larger 
than the action correction in ([75]) . The former one gen- 
erates a shift in the Gaussian integral over j/o which after 
performing the integration leads to a contribution with 
(^„_^g 2/„(7_„)2 ^ l/7'^. The essence of this analysis is 
that if only the first sum in 5a is kept, the distribution 
is calculated with corrections at most of order ln(7)/72. 
This in turn means that the correction 5S in the action 
must be neglected. This way, one again ends up with in- 
dependent Gaussian integrals over the fluctuations modes 
yn, the local frequency A{q) of which reads in this order 
of perturbation theory 



A^V" - pV'V" [Ao - pV'Xl] 



(80) 



The result for the scaled fluctuation factor (PT|) i^(2) 
follows thus from F^^-' in (|53p with the replacement 
V"{q) — > A{q). Contributions to the full path integral be- 
yond the Gaussian approximation for the yn are also neg- 
ligible. Hence, the quantum contribution to the thermal 
ZPf^ /P^g"^^ in the range 7^/3 > 1 



distribution y(2) 
reads 



y(2) 



X e 



1 - pV'Xa + p^ 

-[S'-^'>-hpv]/h 



^y"2 
4 



V'V" 



(81) 



where the action contribution ^^^^ — PhV is the exponent 
in ([54|. This result is exact up to order log(7)2/(72Z). 
Higher order correction terms couple the Fourier modes 
in the minimal action path as well as in the fluctuation 
path integral and are analytically no longer accessible. 
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The above distribution determines the Moyal coefR- 
cients of the QSE to be of the form 



D 



^V'{q) {1 + ^ V"{qY - 2V'{q)V"'{q) 



A 



D 



(2) 



1 



1 



2 pV"{q)\a + /32 [lV"{qY + V'{q)V"'{q)\ A^ ' 

(82) 

Now, also the diffusion coefficient carries information 
about the anharnionicity [cf. ([5^ ] Interestingly, the 
structure of the drift ZJ^^-* is similar to that of d'"-^'' with 
2V'V"' replaced by -2V'V"' . This is a direct conse- 
quence of the perturbation theory for the position dis- 
tribution, which reproduces Pl*^' only with the given set 



of d'^^^ , d''2^ . Apparently, the leading anharmonic cor- 
rections to P/3 appear in the fluctuation prefactor only. 
Eventually, we replace in d{q) the second derivative of 
the potential consistently by A{q) so that 



l-Lv"{qfXl 



(83) 



remains unaltered compared to S-^\ Upon inspection 
of ([STjl and (|82p. we find the third derivative of the po- 
tential to appear in the second order term in Aq. This 
reflects the nonlocal feature of the deep quantum regime 
and shows that Ao indeed plays the role of a semiclassi- 
cal parameter similar to Planck's constant for an isolated 
quantum system. In principle, one must now also extract 
dynamical corrections from the reduced dynamics. How- 
ever, we know from the previous section that in local 
harmonic approximation these are at most of order 1 /^"^ 
compared to the leading term 1/7717. The quantum cor- 
rections accounted for in the above coefficients are much 
larger, namely, of order \og{"f)'^ / (7^/). As a consequence, 
additional dynamical correction due to the anharmonic- 
ity do not play a role. 

To illustrate the perturbative results for the equilib- 
rium distribution, we consider a double well potential 
of the form V{q) = -mujlq^/2 + aq^/4,. In Fig. [3] the 
thermal density is shown in the classical regime together 
with the quantum distributions in leading order, local 
harmonic approximation, and beyond as discussed above. 
The leading order describes the quantum tunneling in the 
barrier range, while the higher order approximations also 
capture the tunneling at the rising walls of the potential. 
Around the barrier top all quantum results basically co- 
incide and start to deviate only towards the well region. 
Further, the maximum of the distribution shifts slightly 
towards the barrier top thus reflecting the finite trans- 
parency of the barrier. 



V. COMPARISON WITH OTHER 
APPROACHES 

The analysis of the strong friction limit in quantum 
mechanics based on path integrals as outlined in 13] has 




X 

FIG. 3: Equilibrium distribution scaled with Pj^^\x = 0) in a 

double well potential vs. the scaled position x — q/^h/mojo 
for various orders of perturbation theory. The classical result 
(dotted) is depicted together with the leading order expres- 
sion (short-dashed) according to (|47l) . the expression in local 
harmonic approximation (dashed) according to (|54|) . and the 
expression beyond (solid) according to H81[) . Parameters are 
7/ajo ~ 3, ujofiP ~ 1 for the bath and ah/m^ujQ — 0.5 for the 
potential so that the well minimum is located at a; = %/2. 



triggered other studies proposing alternative approaches 
to derive quantum Smoluchowski equations. Here, we 
briefly discuss them in comparison with this previous 
work and its extension presented above. For this purpose 
we recall that one carefully has to distinguish between the 
strong friction range at high temperatures 7/1/? <^ 1 and 
that at low temperatures 7?i/3 ^ 1. While in the latter 
domain the strong friction analysis must be based on a 
formulation which captures the system-bath interaction 
non-perturbatively, in the former region one could hope 
that perturbative formulations like e.g. master equations 
are sufficient to derive at least leading quantum correc- 
tions. This assumption is based on the fact that the 
Wigner representation of e.g. the Caldeira-Leggett mas- 
ter equation reduces in the high temperature limit to the 
classical Fokker-Planck equation. 

One approach put forward in [s^j starts directly from 
the classical Smoluchowski equation and quantizes this 
reduced equation of motion. Such a procedure is ques- 
tionable already at sufficiently elevated temperatures and 
it is certainly not a consistent way to take into account 
the quantum mechanics of the reservoir and its inter- 
action with the system. In fact, corresponding results 
contradict the fluctuation-dissipation theorem and may 
thus even lead to unphysical predictions. This failure has 
already been discussed in [3J|. 



The formulation developed in [3J| focuses on the case of 
a QSE for free Brownian motion and uses an ad hoc pro- 
cedure to derive quantum corrections in the mean square 
displacement in position. From the exact path integral 
results (5|,i20|] one finds that in this case {[q{t) — q{0)]'^) = 



13 



2{t + At)/M[3j, where At captures deviations from the 
classical behavior. For strong friction and on the coarse 
grained time scale h/3,l/j <C t [see ()15|) ] At is at most 
of order h/3 and must thus be discarded in a systematic 
treatment. As shown above, the QSE for a free Brow- 
nian particle coincides with the classical Smoluchowski 
equation. 

In the hi gh temp erature domain 7ft/? <C 1 Coffey and 
co-workers [22, [2J| have recently derived quantum cor- 
rections in a QSE by determining a current operator of 
the form ([23|l . however, using the thermal Wigner func- 
tion of the bare system. The crucial question is, whether 
in this domain quantum corrections can then be derived 
systematically. Unfortunately, and in contrast to a naive 
expectation, this is not the case. To see this in detail, 
we look at the exactly solvable harmonic system treated 
already in Sec. IIIII In the strong friction limit the ex- 
act diffusion term in the QSE is given by the variance in 
position TOW§((7^) as shown in (|22|) . i.e.. 
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For jhfi <C 1 this expression can be expanded in powers 
of h. In leading order, only the n = contribution must 
be taken into account and one regains the classical result. 
In next order, only the i^^,n ^ term in the denomina- 
tor must be kept due to <^ 1, which in turn implies 
LOohp <C 1 since 7/wo 3> 1. Hence, the leading quantum 
correction is (uj^ / 13) J2nM'^ / '^n) = t^lh^P/l'i.. This is 
the well-known universal quantum correction, also dis- 
cussed in textbooks which reveals that in the high 
temperature limit the leading quantum correction is in- 
dependent of friction. Hence, in principle it can indeed 
also be derived from the hare equilibrium distribution. 

Now let us look at higher order quantum corrections. 
An ?i-expansion of the Wigner function of the hare system 
as done in [2^ [2^ leads to a powers series exactly with 
(/3/2) X fi^ (3/ 12 as expansion parameter. Accordingly, 
quantum corrections in a diffusion coefficient gained with 
this expansion are independent of 7 to all orders of h. 
This can only be true though if the Vn"f term in the above 
sum can be neglected in all orders in h. This leads to 
the conditions Vnl f^,Wg,n ^ 0. The first relation 
brings us back to ^h(3 1, while the second one gives 
^/ujQ ^ uJoh(3. For strong friction j/ojq ^ 1 these two 
conditions contradict each other since the first one im- 
plies uJoh(3 <^ 1, while the second one requires 1 ^ woft/?. 
It is thus not consistent to neglect the 7-dependence of 
the diffusion term beyond the leading quantum correc- 
tion. The correct expansion in the strong friction range 
which only requires 7/^9 ^ ensure the validity of 

the Markov approximation and applies from ^hp <C 1 to 
7?i/3 ^ 1 is given in ([^5]) . In the high temperature range 
uji^hP <^ jh/S <^ 1 this leads to the series 
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FIG. 4; Quantum corrections SD'^P'^ (|87I) in the diffusion co- 
efficient for a harmonic system for n = 1 (thin) and n = 2 
(thick) . The systematic strong friction expansion according to 
(I28p (solid) is shown together with the "without-7-expansion" 
(23| (dotted). The parameter 7^/3 = 0.4 is kept fixed, while 
7/a;o is varied. For n = 1 both expansions coincide (see text). 



in contrast to Coffey's result 
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The first conclusion is that Coffey's approach is correct 
in the regime 7/1/? <C 1 only in leading order in fi, where it 
coincides with the leading order result presented above. 
Hig her order corrections in the QSE as specified e.g. in 
[23| are incorrect and, while they may give qualitatively 
reasonable results, lead to uncontrolled approximations. 
The second and more fundamental conclusion is that for 
larger friction sub-leading quantum fluctuations are de- 
termined by the system-bath interaction so that any ap- 
proach which treats this interaction pcrturbatively fails 
even at elevated temperatures. 

In the above discussion we focused on the position vari- 
ance for a purely ohmic spectral density. In a strict sense, 
however, the momentum variance diverges in this limit 
[cf. discussion below (pn|) ]. This problem is usually cured 
by introducing a high frequency cut-off Wc as mentioned 
at the end of Sec. Ill Al [5]. The above high temperature 
expansion remains then valid apart from corrections of 
order ^/ujc,(^o/i-^c ff 3> 1. In particular, this lat- 

ter relation guarantees that equilibrium fluctuations are 
still determined by almost ohmic spectral densities and 
almost ohmic friction functions 7(2) in accordance with 
an effectively Markovian theory. 

To illustrate this discussion we show in fig. 2] succes- 
sive quantum corrections in the exact diffusion coefficient 
D^^'' of a harmonic system according to (pS)) and accord- 
ing to Coffey's "without-7-expansion", respectively. It is 
convenient to consider the respective nth order approxi- 
mation I? 2"^ properly weighted by the full result, i.e.. 
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For n = 1 we have = of^ and D^^^ = + 

wlh^P/U so that both expansions coincide as pointed 
out above. Substantial differences appear, however, for 
n = 2 in the range, where the authors of [22l | claim their 
approach to work, namely, for 7/1/3 <C 1 and ^/loq ^ 
1. Even in the strongly underdamped regime for fixed 
7/1/3 ^ 1 the "without-7-expansion" is only of limited 
use since one then leaves the classical regime and enters 
the low temperature range. 

The escape rate in the high temperature range ob- 
tained within this type of approximation is shown in 
fig-Ill According to the above discussion the correspond- 
ing quantum enhancement factor is found to be inde- 
pendent of friction [2^. Obviously, deviations tend to 
increase for stronger friction, i.e. in the regime where the 
approach is supposed to apply. 

VI. DISCUSSION 

In this paper we developed a systematic semiclassical 
expansion for the imaginary time path integral of the 
statistical operator in order to derive drift and diffusion 
terms in the QSE. Since in the QSR the thermal time 
scale hp is still small compared to the time scale for re- 
laxation of the marginal distribution in position, non- 
Markovian features of the reduced dynamics do not play 
a role and a QSE exists. Drift and diffusion are directly 
related to the action of the extremal paths and fluctua- 
tions around them. The physical picture behind is that 
for strong friction and on the thermal time scale hp orbits 
explore the potential only locally. This allows for a per- 
turbative treatment which successively takes higher order 
terms in the expansion of the potential around the aver- 
age minimal action orbit into account. Hence, in leading 



order the dynamics takes effectively place in a potential 
of constant force, in next order in a local harmonic poten- 
tial, and anharmonicities become relevant only beyond. 
Dynamical corrections are classical in nature, at least up 
to the local harmonic approximation. 

The systematically improved QSE is applied to calcu- 
late the escape rate out of a metastable state for low 
temperatures, but above crossover. It then reproduces 
the quantum enhancement factor exactly, whereas dy- 
namical corrections only appear in the classical form of 
the Grote-Hynes frequency. Note that for strong fric- 
tion the inverse crossover temperature /3o is given by 
/i/3o ~ 27r7/a;^. On the other hand the QSR condition re- 
quires hj3 <g; 7/w^. This is why the deep tunneling regime 
below the crossover lies outside the range of validity of 
the QSE and can thus not be captured by this type of 
approximation. Above the crossover, however, the rate 
expression is not sensitive to local anharmonicities of the 
metastable potential, since such contributions only enter 
in third and higher order terms of a variational expansion 
around the trivial minimal action paths residing at the 
well bottom and at the barrier top. 

The systematic expansion outlined here allows to study 
quantum effects in overdamped systems in a broad range 
of parameters comprising the domains of high as well as 
of low temperatures and due to the high order of pertur- 
bation theory also the range of strong to only moderate 
dissipation. 



Acknowledgments 

We thank H. Grabert and A. Verso for fruitful discus- 
sions. 



[1] H. Risken, The Fokker-Planck Equation (Springer, Berlin 

Heidelberg, 1989), 2nd ed. 
[2] R. K. Wangsness and F. Bloch, Phys. Rev. 89, 728 

(1953). 

[3] A. G. Redfield, IBM J. Res. Dev. 1, 19 (1957). 
[4] C. W. Gardiner, Quantum Noise (Springer, Berlin, 
1991). 

[5] U. Weiss, Quantum Dissipative Systems (World Scien- 
tific, Singapore, 2008), 3rd ed. 
[6] H.-P. Breuer and F. Petruccione, The Theory of Open 

Quantum Systems (Clarendon Press, Oxford, 2006). 
[7] V. Averin, B. Ruggiero, and P. Silvestrini, Macroscopic 

Quantum Coherence and Quantum Computing (Kluver 

Academic, Plenum, New York, 2001). 
[8] L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni, 

Rev. Mod. Phys. 70, 223 (1998). 
[9] A. Garg, J. N. Onuchic, and V. Ambegaokar, J. Chem. 

Phys. 83, 4491 (1985). 
[10] J. Jortner and M. B. (eds.). Adv. Chem. Phys. 106, 107 

(1999). 

[11] P. Reimann, Physics Reports 361, 57 (2002). 



[12] M. Smoluchowski, Ann. Phys. 353, 1103 (1915). 

[13] J. Ankerhold, P. Pechukas, and H. Grabert, Phys. Rev. 

Lett. 87, 086802 (2001); ihid 101, 119903 (2008). 
[14] L. Machura, M. Kostur, P. Hanggi, P. Talkner, and 

J. Luczka, Phys. Rev. E 70, 031107 (2004). 
[15] J. Ankerhold and H. Lehle, J. Chem. Phys. 120, 1436 

(2004). 

[16] J. Ankerhold, Europhys. Lett. 67, 280 (2004). 

[17] M.-L. Zhang, S. Zhang, and E. PoUak, J. Chem. Phys. 

120, 9630 (2004). 
[18] L. Machura, M. Kostur, P. Talkner, J. Luczka, and 

P. Hanggi, Phys. Rev. E 73, 031105 (2006). 
[19] J. Dajka, L. Machura, S. Rogozinski, and J. Luczka, Acta 

Physica Polonica B 38, 1737 (2007). 
[20] H. Grabert, P. Schramm, and G.-L. Ingold, Phys. Rep. 

168, 115 (1988). 
[21] P. Pechukas, J. Ankerhold, and H. Grabert, Ann. Phys. 

(Leibzig) 9, 794 (2000). 
[22] W. T. Coffey, Y. P. Kalmykov, S. V. Titov, and B. P. 

Mulligan, J. Phys. A 40, F91 (2007). 
[23] W. T. Coffey, Y. P. Kalmykov, S. V. Titov, and L. Cleary, 



15 



Phys. Rev. E 78, 031114 (2008). 
[24] W. T. Cofley, Y. P. Kalmykov, S. V. Titov, and L. Clcary, 

J. Chem. Phys. 131, 084101 (2009). 
[25] A. O. Caldeira and A. J. Leggett, Physica 121A, 587 

(1983). 

[26] R. Karrlein and H. Grabert, Phys. Rev. E 55, 153 (1997). 
[27] J. L. Skinner and P. G. Wolynes, Physica 96A, 561 
(1979). 

[28] P. Pechukas, J. Ankerhold, and H. Grabert, J. Phys. 
Chem. B 105, 6638 (2001). 



[29] R. Dillenschneider and E. Lutz, Phys. Rev. E 80, 042101 

(2009). 

[30] J. Ankerhold, Quantum Tunneling in Complex Systems 

(Springer, BerUn Heidelberg, 2007). 
[31] P. G. Wolynes, Phys. Rev. Lett. 47, 968 (1981). 
[32] A. O. Bolivar, Phys. Rev. Lett. 94, 026807 (2005). 
[33] J. Ankerhold, H. Grabert, and P. Pechukas, Phys. Rev. 

Lett. 95, 079801 (2005). 
[34] R. Tsekov, J. Phys. A: Math. Theor. 40, 10945 (2007). 



